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We compare the derivation of the dynamic Stark shift of hydrogenic energy levels in a classical 
framework with an adiabatically damped laser-atom interaction, which is equivalent to the Gell- 
Mann-Low-Sucher formula, and a treatment based on time-independent perturbation theory, with 
a second-quantized laser-atom dipole interaction Hamiltonian. Our analysis applies to a laser that 
excites a two-photon transition in atomic hydrogen or in a hydrogenlike ion with low nuclear charge 
number. Our comparisons serve to demonstrate why the dynamic Stark shift may be interpreted 
as a stimulated radiative correction and illustrates connections between the two derivations. The 
simplest of the derivations is the fully quantized approach. The classical and the second-quantized 
treatment are shown to be equivalent in the limit of large photon numbers. 
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I. INTRODUCTION 

The dynamic (AC) Stark shift is a perturbative effect 
that shifts atomic energy levels in a laser field. It is 
an essential topic in precision spectroscopy experiments, 
which have reached unprecedented accuracy^ and gen- 
eral interesti2i^ The dynamic Stark effect is well within 
the reach of students who have studied the quantum me- 
chanics of the hydrogen atom and can be understood at 
different levels in theoretical physics courses. 

The two approaches we will present, the classical field 
and fully quantized field description, coincide in the clas- 
sical limit of high photon density, as is to be expected. 
So the AC Stark effect also serves as an example of how 
to connect and contrast classical and quantum notions of 
physical phenomena. 

If an atom is exposed to external electromagnetic 
fields, its energy levels are shifted due to the interaction 
of the electrons with the field. This shift of energy lev- 
els can be observed in spectroscopic experiments, where 
for example, absorption or fluorescence spectra are mea- 
sured. The Zeeman effect describes this energy shift for 
static magnetic fields, and the DC Stark effect is respon- 
sible for the level shift in static electric fields. Both these 
effects can be avoided in principle by a proper shielding 
of the atom. However, the probing laser light with which 
atoms are irradiated in order to obtain a spectrum also 
constitutes a time-dependent electromagnetic field and 
is necessarily always present in laser spectroscopy. Its 
impact on atomic energy levels is called the AC Stark 
effect and for nondegenerate states it can be understood 
as a time averaged DC Stark shift, as explained in the 
Appendix. This statement holds only for off-resonant 
driving of the atom, where the AC Stark shift can be 
considered as a perturbation. 

In this article off-resonant driving is to be understood 
with regard to any electric-dipole allowed, one-photon 
transition. Even in the nonresonant case two-photon 
transitions can be driven effectively when the frequency 
of the incident radiation is close to half the atomic tran- 
sition frequency. 



Before we discuss the off-resonant excitation of an 
atom by laser radiation, we briefly mention the main dif- 
ferences with the resonant casei^ Most importantly, if 
the frequency of the incident radiation is close to a one- 
photon resonance, then the influence of the laser field 
on the atomic levels lies outside the regime of perturba- 
tion theory and must be included nonperturbatively in 
the dressed-state picture. The reason is that the dipole 
matrix element between the states involved is nonzero 
for one-photon transitions. Consequently, the level shift 
is linear in the electric field amplitude of the laser, in 
contrast to the quadratic dependence that we will obtain 
for off-resonant excitation. For further information on 
dressed states we refer the reader to Rcfs. "s^^ MTollTTl 

In a classical framework the dynamic Stark shift can 
be described by time-dependent perturbation theory. We 
will demonstrate that the dynamic Stark shift can be 
used to illustrate some basic aspects of quantum electro- 
dynamics (QED). The dominant shift of the energy levels 
in this case can be attributed to a second-order pertur- 
bation in which a laser-photon is created or annihilated 
in a virtual intermediate state. 

This article is also devoted to showing that the AC 
Stark shift can be identified as a stimulated radiative 
correction. Indeed, the AC Stark shift is approximately 
equivalent to a spectral component of the electron self- 
energy (Lamb shift) that results when we restrict the dis- 
cussion to a single mode of the electromagnetic field, but 
with an important difference: for the Lamb shift the pho- 
ton modes are all unoccupied in the unperturbed state 
in contrast to the AC Stark shift for which there is one 
highly occupied mode of the electromagnetic field, the 
laser mode. 

This paper is organized as follows: In Sec. the dy- 
namic Stark shift is described using a classical laser field, 
which necessitates the use of time-dependent perturba- 
tion theory and an adiabatically damped interaction. In 
contrast in Sec. IIIII we derive the dynamic Stark shift 
using a quantized-field approach and time-independent 
perturbation theory. The classical and second-quantized 
results are shown to agree in the classical limit, that is, 
for a macroscopically populated laser field mode. 
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II. CLASSICAL FIELD APPROACH 

In this section we rederive the classical expressions for 
the dynamic (AC) Stark shift using a classical description 
of the laser field. Our approach is the usual one employed 
in the literature, and our treatment and our notation 
is inspired by Chap. 5 of Ref. ^3 Let us consider the 
Hamiltonian 



H = Ho + V{e,t), 



where 



Ha = 



2toc Aneor 
V{e,t) = Fexp(-e|i|)cos(tJLt), 
V = -ezSi,. 



(1) 



(2a) 

(2b) 
(2c) 



The Hamiltonian in Eq. describes a hydrogen atom 
(Z = 1) or a hydrogcn-likc ion of nuclear charge number 
Z > 1 in a plane-wave monochromatic laser field, polar- 
ized along the z-direction, adiabatically damped in the 
distant past {t —> —oo) and the distant future {t — > oo). 
Hq describes the nonrelativistic unperturbed hydrogen 
Hamiltonian, and V{e,t) is the time-dependent, adiabat- 
ically damped, harmonic perturbation with magnitude 
V; e is the infinitesimal damping parameter (see, for ex- 
ample, p. 342 of Ref. ilS) . We have assumed that the 
wavelength of the driving light of angular frequency wl 
is large compared to the spatial extent of the atomic 
wave functions (the dipole approximation). The laser- 
atom interaction V{e, t) is treated in the length gauge 
as in Ref. 14 with electric field amplitude £i^. The elec- 
tric field strength involved in V{e, t) is a gauge-invariant 
quantityii^ 

The parameter e > is introduced to avoid a sudden 
turn-on of the perturbation. In the limit e — *■ we will 
obtain the constant intensity result after carrying out the 
relevant time integrations of the first few terms in the 
Dyson series. The introduction of an adiabatic damping 
parameter is also a key element of time-dependent pertur- 
bation theory in QEDi^ In QED the interaction Hamil- 
tonian is usually expressed in the interaction picture and 
a time-dependence is incurred for the field operators (see 
Appendix A of Ref. 6) . 

Energy shifts in QED are usually formulated using the 
Gell-Mann-Low-Sucher theoremiiSiii The applicability of 
this theorem is not restricted to the case of perturba- 
tions in a second-quantized approach, but can be applied 
equally well to a time-dependent, classical perturbation. 

We now consider the effect of the off-resonant pertur- 
bation by an time-dependent electric field on a reference 
state \4>) of the unperturbed atom. In the interaction pic- 
ture (denoted by the subscript /), V{e,t) is represented 
by 

Vi{e,t) = exp (^i/oi)^(e,t)exp ( - '-H^t). (3) 



From the Dyson series we can calculate the time evolution 
operator Ui{e,t) up to second order in Vi: 



C/i(e,t) 



^-Jl I dt'Vi{e,t') 



dt' / dt"Vi{e,t')Vi{e,t").{A) 



Now consider the time-dependent atomic state |'i/'i(^)} hi 
the interaction picture subject to the initial condition 
\ipi{t = —oo)) = \(j)), where the reference state |0} is an 
eigenstate of the unperturbed Hamiltonian Hq. We ex- 
pand IV-'i(O) i'^ ^ complete set {|m)} of eigenstates of Hq 
as 

IV^i(t)) = [/i(e,<)|^i(-oo)) =^c„(t)|TO), (5) 

where Cm{t) — {m\ipi{t)). The initial condition is thus 
C0(— oo) = 1 for the reference state |0} with all other 
c„i(— oo) equal to zero. We are interested in the projec- 
tion 



c4t) = (0|Vi(i))-(0|C/i(e,t)|<^). 



(6) 



We substitute Ui{e,t) from Eq. Q and because {(j)\z\(f)) 
vanishes for parity eigenstates the leading order is 
and the problem reduces to calculating the matrix 
element 



M 



t t' 

dt' [ dt"{(f>\Vi{€,t')Vi{e,t")\(l)) 



t' 



(7a) 



dt' dt"{(l)\Vi{e,t')\m){m\Vi{e,t")\^), 



(7b) 

where the multi-index m counts all bound and contin- 
uum states of the unperturbed hydrogen atom. Because 
the perturbation is harmonic, the time integrals can be 
done without difficulty, convergence being ensured by the 
adiabatic damping. We obtain 



M 



hl^ (0|l/|TO)(m|y|(/.)exp(2ei) 
7 4 ^ 2e(£'0 - E„, ± fiiJL + «fte) ' 



(8) 



with V as defined in Eq. if^ : represents the energy 
of the unperturbed atomic state \(j)). The ± index de- 
notes the summation of the two terms differing only in 
the sign of ftcjL in the denominator. This sum and the 
factor of 1 /4 originate from the definition of the cosine in 
terms of exponential functions. In view of Eqs. Q, lO, 
and © , we have in second-order time-dependent pertur- 
bation theory 



m,± 



(0|y|TO)(m|y|0)exp(2et) 
2e(£'0 - Em ± /iwl + ifie) 



(9) 
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where higher-order terms have been neglected. Now con- 
sider 



d_ » ^ {cl,\V\m){m\V\<l>) 

■m,± ^ 



(10) 



Here the logarithm has been expanded up to second or- 
der in V and exp(2ei) has been replaced by unity. The 
solution of Eq. (|10|l implies that 



c^(t) =cxp(- -^ASAcWi), 



(11) 



where we have defined the dynamic Stark shift A_Eac (0) 
of the reference state 16) 



AEacW = 7^ 



\V\m){m\V\ 



(12) 



In view of Eq. jSJ, we have 

IV'i(i)) =c^WI0) + ---, (13) 

where the ellipsis denotes the projections onto the non- 
reference atomic states. Because the Schrodinger picture 
wave function is related to its interaction-picture coun- 
terpart via |'0(i)} — exp(— ii/oi)|V'i(^)}: we have 



(14) 



The nonreference states from Eq. (|13|l give no contribu- 
tion because they are orthogonal to \4)). The projection 
f(T^ yields the influence of the perturbation on the refer- 
ence state 10) by projecting the time-evolved perturbed 
state onto the reference state such that the perturbation 
to the eigenenergy E^ can be directly seen. 

Note that AEac{4') can in general be complex, rather 
than real. We define 



-f^^--Im{AEAcm, 
AE^^Re{AEAc{<l>))- 



(15) 
(16) 



The real part of the AC Stark effect describes the energy 
shift of the unperturbed energy E^, and the imaginary 
part, if present, can be interpreted as the ionization rate 
70. We can now express the dynamic Stark shift of |0) 
as 



AEacW = \Y. 



V- 



1 



Ed, — Ho ± fkJi^ + ihe 



-V 



(17) 

where the closure relation for the spectrum is employed. 
In the Appendix the zero-frequency limit of Eq. (|17|) is 
related to the static Stark effect. Equation Ijl?!) can be 
written conveniently as a product of a prefactor and a 
sum of two matrix elements, where E is the energy of 
the respective intermediate state E = E^ =p fiioi^: 



AEacW 



± 



(j)\z 



1 



Ho - E^± hwi. 



(18a) 



2ceo 



where {'f) is the dynamic polarizability of the atom in 
the reference state for angular frequency wl of the driving 
laser field. The intensity / of a plane electromagnetic 
wave is 



/ = -eocSl. 



(19) 



This derivation completes our analysis of the AC Stark 
shift using an adiabatically dampedji^ classical-fieldi^ 
approach. 

To illustrate the connection to the Gell-Mann-Low- 
Sucher theorem, we observe that C0(O) in Eqs. JHl and 
@ can be identified with the quantity {a\U{Q, — oo; e)|Q;) 
in the notation of Ref. 0, and the perturbation V{e,t) 
as defined in Eq. (|2bp has to be supplemented by an aux- 
iliary scaling variable g, 



V{e,t) ^ gV{e,t). 



(20) 



The parameter g is later set equal to unity. We then have 
according to Eq. (2±) of Ref.llTL 



AE = lim (ihe g^\n{4>\Ui{€M(l))] (21a) 
g~*i \ ag J 



lim (ihe q-^ In (l — g 



(21b) 



hm^ ^ ^ V (0|gy|m)(m|g^|0) 
g-i 8dg\^Es~ E„, ± hwi^ + ihe ' ^ ' 



h\V\'m)l^m\V\(\>) 
4 ^ E^ - Em ± ftwL + ihe 



(21d) 



The latter result AE = AEac{<I>) agrees with Eq. l(T^ . 
In the step leading to Eq. (|21c|) , an expansion of the log- 
arithm in powers of g is implied, which is equivalent to a 
second-order expansion in the time-dependent perturba- 
tion V. 



III. FULLY QUANTIZED APPROACH 

In the classical picture we set the field amplitude to 
a constant value £l and used time-dependent perturba- 
tion theory with an adiabatic damping parameter. When 
treating the light as a photon field, the classical picture 
can be interpreted as the limit of the fully quantized 
treatment in the limit of a large photon number. 

In second quantization the Hamiltonian for the coupled 
system, atom + radiation field, reads 



(22) 



The first term contains a sum over the discrete spec- 
trum and an integral over the continuous spectrum of 
the Schrodinger equation. We do not consider electron- 
positron pair creation, and therefore we do not quantize 
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the fermion field. The laser field is described as a quan- 
tized photon field with creation and annihilation oper- 
ators and cl, respectively. Hi, reads (in the length 
gauge) 



(OL 



(23) 



(See also Eqs. (4.7) and (4.8) of Ref.i6|.) The symbol V 
denotes the normalization volume and is chosen so that 
the energy density of a one-photon Fock state when in- 
tegrated over V yields fiwL- It might be argued that 
a coherent state of the photon field is a much better 
description than a Fock state with til photons in the 
laser mode, which we have assumed here. However, in 
the limit of large photon number, the relative fluctuation 
of the photon number Sni,/ni, goes to zero for a coher- 
ent state, and we may therefore resort to the Fock-state 
approximation^ 

We work in the Schrodinger picture where the field 
operators carry no time dependence. It is not so widely 
known that it is possible to formulate time-independent 
operators for the quantized radiation field, let alone to do 
meaningful calculations with these operators. However, 
this formulation is introduced in a few textbooks such as 
Ref.iia 

The concept of time-independent field operators has 
also been used for quantum electrodynamic calculations 
(see for example, Eq. (5) of Ref. 1201) . Following this 
approach, we are now in a position to apply time- 
independent perturbation theory^ This approach leads 
to the following second-order result for the energy shift 
of the unperturbed eigenstate |0, n^). 



,ni,\Hi,\m,ni, - l)(m, til - l|gL|0, "-l) 

+ ni^huJi, - {E,n + (n-L - 1)^l) 
(0,nL|i?L|"^, + l)(TO,nL + l\Hi,\(j), ni,) 



E(j, + TlLfiaJL ~ (Em + [riL + l)fi^L) 
{(j)\z\m) (to|z| 



2eoV lEcf, — E„ 
{4>\z\m){m\z\ 



huji 



Ea, — E„ 



(^L + 1) 



(24) 



The sum over virtual intermediate states \m) has been 
split into two parts depending on the number of photons 
in the field. In the classical limit riL — > cxd, V — > oo, 
riL/V = constant, we can simplify Eq. (|24|) to read 



A£^Ac(0) 



?^JmiML 
2eoV 



(25) 



with Pi^L (0) given in Eq. H18a|) . 

The remaining issue concerns the matching of this re- 
sult to the classical result in Eq. p8|l . In the quantized 
formalism the term 



w 



ni^hioi 
V 



(26) 



gives the energy density in which the atom is immersed, 
which is related to the intensity via 



/ = wc. 

We use Eqs. and obtain 



Ai;Ac(<^) 



2enc 



(27) 



(28) 



in agreement with Eq. 1)18(1 . Thus the classical-field and 
the quantized-field approach give consistent results in the 
classical limit. 



IV. CONCLUSIONS 

We have contrasted two ways of deriving analytic ex- 
pressions for the dynamic Stark shift of a hydrogenic en- 
ergy level. The first, based on an adiabatically damped 
length-gauge interaction (see Eq. (|2bl) '). leads to a clas- 
sical treatment where the electric laser field is simply 
modeled as a periodic perturbation (see Sec. |n}. The 
second derivation, based on a quantized description of 
the electromagnetic field, leads to expressions that are 
equivalent to the classical expressions in the limit of a 
large occupation number of the laser mode (see Sec. lIIIIl . 

The AC Stark shift has been characterized as a stim- 
ulated radiative correction—- because it results from a 
self-energy-like formalism if the sum over virtual modes 
of the photon field is restricted to a single mode: the 
laser mode. We illustrated this statement by giving an 
explicit derivation in Eq. I|24|l . This treatment is based 
on time-independent field operators. Equation l|24(l illus- 
trates how the classical predictions should be modified in 
an environment where the photon number is not large. 
Indeed, the AC Stark shift receives an interpretation in 
this context as the second-order perturbation incurred 
by the coupled system, atom -f laser field, due to vir- 
tual creation and annihilation of laser photons. When 
the perturbation is evaluated using an empty Fock space 
as the unperturbed state, and when a sum is formed 
over all possible virtual excitations, the self-energy is 
obtainediS*2i. 
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APPENDIX A: DC STARK SHIFT 



We use the relation 



For nondegenerate states, the DC Stark shift is a 
second-order perturbation in the electric field strength 
and can be interpreted as the zero-frequency limit of the 
AC Stark shift. In this appendix we briefly illustrate 
this relation. In nonrelativistic quantum mechanics the 
ground state is the only nondegenerate state. However, 
as a consequence of the spin of the electron (fine struc- 
ture), the spin of the nucleus (hyperfine structure) and 
QED effects (Lamb shift), the degeneracy of other states 
is broken, and as long as the DC Stark shift is small 
compared to these energy differences, the following con- 
siderations also hold for excited states. For larger per- 
turbations, the DC Stark effect is linear in the electric 
field. 

For a state |0) that fulfills the above conditions, con- 
sider the limit of the dynamic Stark shift obtained in 
Eq. (|I8|I as the angular frequency of the laser field goes 
to zero 



lim — 

ri^Q 2 

-t 



1 



\m)(m\ 



77 Ho - E^-ri 



1 



Ho — Ed 



Er, 



Ea, 



(A2) 



where the reduced Green function HA2|I excludes the ref- 
erence state for which the denominator would diverge to 
obtain 



/^E 



AC,0 



fc- C-l 



\ Ho — Ed 



AE- 



DCi 



(A3) 

which is the expression for the second order DC Stark 
shift. The static electric field strength ^dc is matched to 
the harmonic laser field EL{t) — cos(wLi) by averaging 
the laser field strength squared over one optical period: 



AE. 



ACfi 



lim 



4 A. 



(f)\z 



1 



Ho - Ed, ± hwv 



(Al) 



1 



Elit) = -£i 



-DC- 



(A4) 
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